The Physics of Wind-Fed Accretion 



Christopher W. Mauche*, Duane A. LiedahF, Shizuka Akiyama + and 

Tomasz Plewa** 

* Lawrence Livermore National Laboratory, L-473, 7000 East Avenue, Livermore, CA 94550 
^KIPAC, Stanford University, 2575 Sand Hill Road, M/S 29, Menlo Park, CA 94025 
** Florida State University, Department of Scientific Computing, DSL 443, Tallahassee, FL 32306 

Abstract. 

We provide a brief review of the physical processes behind the radiative driving of the winds 
of OB stars and the Bondi-Hoyle-Lyttleton capture and accretion of a fraction of the stellar wind 
by a compact object, typically a neutron star, in detached high-mass X-ray binaries (HMXBs). In 
addition, we describe a program to develop global models of the radiatively-driven photoionized 
winds and accretion flows of HMXBs, with particular attention to the prototypical system Vela 
X-l. The models combine XSTAR photoionization calculations, HULLAC emission models appro- 
priate to X-ray photoionized plasmas, improved models of the radiative driving of photoionized 
winds, FLASH time-dependent adaptive-mesh hydrodynamics calculations, and Monte Carlo radi- 
ation transport. We present two- and three-dimensional maps of the density, temperature, velocity, 
ionization parameter, and emissivity distributions of representative X-ray emission lines, as well 
as synthetic global Monte Carlo X-ray spectra. Such models help to better constrain the properties 
of the winds of HMXBs, which bear on such fundamental questions as the long-term evolution of 
these binaries and the chemical enrichment of the interstellar medium. 
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INTRODUCTION 

As described by Castor, Abbott, and Klein (hereafter CAK) [1], mass loss in the form 
of a high- velocity wind is driven from the surface of OB stars by radiation pressure on a 
multitude of resonance transitions of intermediate charge states of cosmically abundant 
elements. The wind is characterized by a mass-loss rate M ps 10 -10 -5 M Q yr _1 and a 
velocity profile v(r) ps Voo(1 —R*/r)P, where r is the distance from the OB star, /3 ps i, 

Voo ~ 3v esc is the terminal velocity, v esc = (2GM+/R*) 1 / 2 is the escape velocity, and 
and R± are respectively the mass and radius of the OB star. As described by Hoyle and 
Lyttleton [2], Bondi [3], and Shapiro and Lightman [4], in detached high-mass X-ray 
binaries (HMXBs), a compact object, typically a neutron star, captures the fraction of 
the wind of the OB star that passes within the accretion radius r a = £2GM x /v 2 el of the 

compact object, where £ is a constant of order unity, v re i = (v^ + the velocity 

of the wind relative to the compact object, v w is the velocity of the wind at r = a, v x = 
Ina/Po^ is the orbital velocity of the compact object, a = [G(M^ + M x )(i , orb /2^:) 2 ] 1 / 3 
is the binary separation, and P ^ is the binary orbital period. The mass-accretion rate 
onto the compact object M a = 7rr 2 p w v re i, where p w = M w /4na 2 v w is the wind mass 
density and M w is the wind mass-loss rate, hence M a = \ (r a /a) 2 (v re i/v w )M w = fM w . 
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This compact expression expands out to 

M a » O.O12C 2 T]-^M x /M ) 2 (M^M )- 8 / 3 (^/ J R ) 2 ( J P orb /d)- 4 / 3 M w , (1) 

where r) = (v w /v esc ) ~ 3 (1 —R*/a)P for an isolated OB star. For the parameters ap- 
propriate to the prototypical detached HMXB Vela X-l (§2), ! M a « 7 x 10~ 5 M W ~ 
7 x 10~ H M yr . Accretion of this material onto the neutron star powers an X-ray 
luminosity L x = GM SL M X /R X ~ 1 x 10 36 erg s -1 , where M x and 7? x are respectively the 
mass and radius of the compact object. The resulting X-ray flux photoionizes the wind 
and reduces its ability to be radiatively driven, both because the higher ionization state of 
the plasma results in a reduction in the number of resonance transitions, and because the 
energy of the transitions shifts to shorter wavelengths where the overlap with the stellar 
continuum is lower. To first order, the lower radiative driving results in a reduced wind 
velocity v w near the compact object, which increases the accretion radius r a , which in- 
creases the accretion efficiency /, which increases the X-ray luminosity L x . In this way, 
the X-ray emission of HMXBs is the result of a complex interplay between the radiative 
driving of the wind of the OB star and the photoionization of the wind by the neutron 
star. 

Known since the early days of X-ray astronomy, HMXBs have been extensively stud- 
ied observationally, theoretically [5, 6, 7], and computationally [8, 9, 10, 11]. They are 
excellent targets for X-ray spectroscopic observations because the large covering frac- 
tion of the wind and the moderate X-ray luminosities result in large volumes of photoion- 
ized plasma that produce strong recombination lines and narrow radiative recombination 
continua of H- and He-like ions, as well as fluorescent lines from lower charge states. 

VELA X-l 

The prototypical detached HMXB Vela X-l has been studied extensively in nearly 
every waveband, particularly in X-rays, since its discovery as an X-ray source during 
a rocket flight four decades ago. It consists of a BO. 5 lb supergiant and a magnetic 
neutron star in an 8.964-day orbit. From an X-ray spectroscopic point of view, Vela X- 
1 distinguished itself in 1994 when Nagase et al. [12], using ASCA SIS data, showed 
that, in addition to the well-known 6.4 keV Fe Ka emission line, the eclipse X-ray 
spectrum is dominated by recombination lines and continua of H- and He-like Ne, 
Mg, Si, S, Ar, and Fe. These data were subsequently modeled in detail by Sako et al. 
[13], using a kinematic model in which the photoionized wind was characterized by the 
ionization parameter t, = L x /nR 2 , where R is the distance from the neutron star and n 
is the number density, given by the mass-loss rate and velocity law of an undisturbed 
CAK wind [specifically, n = M w /4niiv(r)r 2 , where /i is the mean atomic weight and 
v( r ) = vo + Voo(l —R*/r)P]. Vela X-l was subsequently observed with the Chandra High 
Energy Transmission Grating (HETG) in 2000 for 30 ks in eclipse [14] and in 2001 for 



1 M+ = 23.8 M Q , = 30 R®, M x = 1.86 M , P orb = 8.964 d, Voc = 1700 km s"\ J3 = 0.8, and 
M w w 10~ 6 M @ yr"\ hence a = 3.7 x 10 12 cm, r a = 5.7 x 10 10 cm, v x = 300 km s" 1 , v w = 880 km s"\ 
vw — 930 km s _1 , v esc = 550 km s _1 , and 77 = 1.6. 
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FIGURE 1. Three views of the Chandra HETG Medium Energy Grating count spectrum of Vela X- 1 in 
eclipse, with the strongest emission features labeled: (a) the broad 0.7-10 keV bandpass, (b) the line-rich 
4-15 A bandpass, and (c) the 6.1-7.3 A bandpass containing emission lines from H- and He-like Si as 
well as fluorescence emission features from L-shell Si. 



85, 30, and 30 ks in eclipse and at binary phases 0.25 and 0.5, respectively [15, 16]; three 
views of the line-rich eclipse spectrum are shown in Figure 1. Watanabe et al. [16], using 
assumptions very similar to those of Sako et al. and a Monte Carlo radiation transfer 
code, produced a global model of Vela X-l that simultaneously fits the HETG spectra 
from the three binary phases with a wind mass-loss rate M w »2x 10~ 6 M Q yr _1 and 
terminal velocity Vc*. = 1 100 km s . One of the failures of this model was the velocity 
shifts of the emission lines between eclipse and binary phase 0.5, which were observed to 
be Av 400-500 km s _1 , while the model simulations predicted Av ~ 1000 km s _1 . In 
order to resolve this discrepancy, Watanabe et al. performed a 1 -dimensional calculation 
to estimate the wind velocity profile along the line of centers between the two stars, 
accounting, in an approximate way, for the reduction of the radiative driving due to 
photoionization. They found that the velocity of the wind near the neutron star is lower 
by a factor of 2-3 relative to that of an undisturbed CAK wind, which was sufficient 
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to explain the observations. However, these results were not fed back into their global 
model to determine the effect on the X-ray spectra. Equation 1 suggests that, for a given 
wind mass-loss rate, the X-ray luminosity should be higher by a factor of (2-3) 4 = 16- 
81, or, for a given X-ray luminosity, that the wind mass-loss rate should be lower by a 
similar factor. 

HYDRODYNAMIC SIMULATIONS 

To make additional progress in our understanding of the wind and accretion flow of 
HMXBs in general and Vela X-l in particular — to combine the best features of the 
hydrodynamic models of Blondin et al. [8, 9, 10, 1 1] and the kinematic- spectral models 
of Sako et al. [13] and Watanabe et al. [16] — we have undertaken a project to develop 
improved models of radiatively-driven photoionized accretion flows, with the goal of 
producing synthetic X-ray spectra that possess a level of detail commensurate with the 
grating spectra returned by Chandra and XMM-Newton. This project combines: 

• XSTAR [17] photoionization calculations, which provide the heating T and cooling 
A rates of the plasma and the ionization fractions of the various ions as a function 
of temperature T and ionization parameter t, . 

• HULLAC [18] emission models appropriate to X-ray photoionized plasmas. 

• Improved models of the radiative driving of the photoionized wind. 

• FLASH [19] two- and three-dimensional time-dependent adaptive-mesh hydrody- 
namics calculations. 

• Monte Carlo radiation transport [20]. 

Radiative driving of the wind is accounted for via the force multiplier formalism [1], 
wherein the radiative acceleration g md = g e (1 + M), where g e = KzL*/4nr 2 c is the ra- 
diative acceleration due to electron (Thomson) scattering; fQ = n e Oj/p fh 0.34 cm 2 g _1 
is the Thomson opacity; n e is the electron number density; p is the mass density; 
M = Iiines( i7 v/i 7 )(v t hV/c)(l - e-W)/t is the line force multiplier; v th = {2kT*/m v ) l l 2 
is the mean thermal velocity of the protons; T+, L+, and F v /F are respectively the 
effective temperature, luminosity, and normalized flux distribution of the OB star; 
t = K e v t hp(dv/dr)~ l is the optical depth parameter; dv/dr is the velocity gradient; r\ 
accounts for all the atomic physics (wavelengths, oscillator strengths, statistical weights, 
and fractional level populations); and the other symbols have their usual meanings. We 
calculated M[T, £ ,/] on a [2 1x51x41] lattice accounting for X-ray photoionization and 
non-LTE population kinetics using HULLAC atomic data for 2 x 10 6 lines of 35,000 
energy levels of 166 ions of the 13 most cosmically abundant elements. Compared to 
Stevens and Kallman [7], our line force multipliers are typically higher at low t because 
of the larger number of lines in our models, and lower at high t because we do not 
assume that the level populations are in LTE. 

In addition to the usual hydrodynamic quantities, the FLASH calculations account 
for: 

• The gravity of the OB star and neutron star (although, following Ruffert [21], for 
numerical reasons we softened the gravitational potential of the neutron star over a 
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FIGURE 2. Color-coded maps of (a) logT(K) = [4.4,8.3], (b) log«(cm- 3 ) = [7.4,10.8], (c) 
logv(kms _1 ) = [1.3,3.5], and (d) log<i; (ergcms -1 ) = [1.1,7.7] in the orbital plane of Vela X-l. The 
positions of the B star and neutron star are shown by the circle and the "x," respectively. The horizontal 
axis x = [—5,7] x 10 12 cm and the vertical axis y = [—4,8] x 10 12 cm. 

scale length e = 10 3 i? x = 10 9 cm). 

• Coriolis and centrifugal forces. 

• Radiative driving of the wind as a function of the local ionization parameter, 
temperature, and optical depth. 

• Photoionization and Compton heating of the irradiated wind. 

• Radiative cooling of the irradiated wind and the "shadow wind" behind the OB star. 

2D HYDRODYNAMIC SIMULATIONS 

To demonstrate typical results of our simulations, we show in Fig. 2 color-coded maps of 
the log of the (a) temperature T, (b) density n, (c) velocity v, and (d) ionization parameter 
^ of a FLASH 2-dimensional simulation in the binary orbital plane of an HMXB with 
parameters appropriate to Vela X-l (specifically, those of Sako et al. [13]). In this 
simulation, the computational volume x = [—5, 7] x 10 12 cm and y = [—4, 8] x 10 12 cm, 
the B star is centered at [x,y] = [0, 0], the neutron star is centered at [x,y] = [0, 3.7 x 10 12 ] 
cm, and 128 x 128 = 16,384 computational cells were used for a spatial resolution 
of Al = 9.4 x 10 10 cm. At the time step shown (t = 100 ks), the relatively slow (v « 
400 km s -1 ) 2 irradiated wind has reached just ~ 2 stellar radii from the stellar surface. 
The various panels show (1) the effect of the Coriolis and centrifugal forces, which cause 
the flow to curve clockwise, (2) the cool, fast wind behind the B star, (3) the hot, slow, 
irradiated wind, (4) the hot, low-density, high-velocity flow downstream of the neutron 
star, and (5) the bow shock and two flanking shocks formed where the irradiated wind 
collides with the hot disturbed flow in front and downstream of the neutron star. 

Given these maps, it is straightforward to determine where in the binary the X- 
ray emission originates. To demonstrate this, we show in Fig. 3 color-coded maps of 
the log of the emissivity of (a) SiXIV Lya, (b) SiXIII Hea, (c) FeXXVI Lya, and 



2 Note that this velocity reproduces the value that Watanabe et al. found was needed to match the velocity 
of the emission lines in the Chandra HETG spectra of Vela X-l. 
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FIGURE 3. Similar to Fig. 2, but for the log of the X-ray emissivity of (a) SiXIV Lya, (b) SiXIII 
Hea, (c) FeXXVI Lya, and (d) FeXXV Hea in the orbital plane of Vela X-l. In each case, two orders 
of magnitude are plotted. 



(d) FeXXV Hea. The gross properties of these maps are consistent with Fig. 24 of 
Watanabe et al., but they are now (1) quantitative rather than qualitative and (2) specific 
to individual transitions of individual ions. The maps also capture features that otherwise 
would not have been supposed, such as the excess emission in the H- and He-like Si 
lines downstream of the flanking shocks. Combining these maps with the velocity map 
(Fig. 2c), these models make very specific predictions about the intensity of the emission 
features, where they originate, and their velocity widths and amplitudes as a function of 
binary phase. 

3D HYDRODYNAMIC SIMULATIONS 

Having developed and extensively tested our software on various 2-dimensional sim- 
ulations, we next calculated a more limited number of 3-dimensional simulations of 
an HMXB with parameters appropriate to Vela X-l. Figure 4 shows, for one of these 
simulations, the log of the number density on three orthogonal planes passing through 
the computational volume, which spans Ax = 12 x 10 12 cm, Ay = 16 x 10 12 cm, and 
Az = 8 x 10 12 cm using 192 x 256 x 128 = 6.3 x 10 6 computational cells, for a spa- 
tial resolution Al = 6.3 x 10 10 cm. In order to compare our hydrodynamic calculations 
to the kinematic models of Sako et al. [13] and Watanabe et al. [16], we also calcu- 
lated, on the same 3-dimensional grid, using the parameters of Watanabe et al., the ve- 
locity v{r) = vq + v<>o(l —R*/r)P , density n = A/ w /47T/iv(r)r 2 , and ionization parameter 
% =L x /nR 2 for an assumed undisturbed CAK wind, assuming that the temperature T ) 
is that of a photoionized plasma in thermal equilibrium: T(T,£) = A(T,£). Figure 1 of 
Akiyama et al. [22] shows color-coded maps in the binary orbital plane of the density 
and temperature and the SiXIV and FeXXVI ion abundances for the CAK and FLASH 
models, emphasizing the dramatic differences between them. 

SYNTHETIC X-RAY SPECTRA 

In the next step in this process, we fed the physical parameters of the CAK and FLASH 
models into our Monte Carlo radiation transfer code and followed the spatial and spectral 
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FIGURE 4. Color-coded maps of logn (cm -3 ) = [7.4, 12. 1] on three orthogonal planes from a FLASH 
3D simulation of Vela X-l, showing the B star (blue sphere), neutron star, and accretion wake. The 
computational volume x = [—6,6] x 10 12 cm, y— [—6,10] x 10 12 cm, andz= [-4,4] x 10 12 cm. 



evolution of photons launched from the surface of the neutron star until they are either 
destroyed or escape the computational volume. As detailed by Mauche et al. [20], the 
Monte Carlo code accounts for Compton and photoelectric opacity of 446 subshells of 
140 ions of the 12 most abundant elements. Following photoabsorption by K-shell ions, 
we generate radiative recombination continua (RRC) and recombination line cascades 
in a probabilistic manner using the recombination cascade calculations described by 
Sako et al. [13], with the shapes of the RRCs determined by the functional form of the 
photoionization cross sections and the local electron temperature. Fluorescence emission 
from L-shell ions is ignored, since it is assumed to be suppressed by resonant Auger 
destruction (however, see Fig. 1(c) and Liedahl [23]). 

The synthetic spectra resulting from the Monte Carlo simulations are shown in Fig. 5. 
The upper histogram is the assumed (featureless) power-law spectrum emitted by the 
neutron star, while the synthetic spectra for the CAK and FLASH physical models are 
shown by the middle and lower histograms, respectively. While it is too early in the 
development of our program to draw any firm conclusions from these spectra (or from 
the differences between them), they demonstrate the type of results produced by our 
Monte Carlo code — X-ray spectra rich in emission lines and RRCs — and the type of 
differences between simple CAK and detailed hydrodynamic models, highlighting the 
importance of the details of the 3-dimensional model to the emitted spectra. 

We close by noting that Fig. 5 plots global X-ray spectra, irrespective of viewing 
angle or occultation by the B star, although it is straightforward to window the output of 
the Monte Carlo code to produce synthetic spectra for a specific viewing direction (for a 
given binary inclination and binary phase) or for a given binary inclination and a range 
of binary phases, and hence allow direct comparisons to actual data. The next step in our 
program is a careful comparison of such synthetic spectra to the grating spectra of Vela 
X-l and other HMXBs returned by Chandra and XMM-Newton. In this way, it will be 
possible to better constrain the various parameters of the winds of HMXBs, such as the 
mass-loss rate M w , terminal velocity v^, velocity profile v(r), and elemental abundances 
of the OB star — parameters that bear on such fundamental questions as the long-term 
evolution of these binaries and the chemical enrichment of the interstellar medium. 
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FIGURE 5. Monte Carlo X-ray spectra of the CAK (red) and FLASH (blue) models of Vela X-l for 
the assumed power law neutron star X-ray spectrum (black). The H- and He-like ions responsible for the 
strongest emission lines and radiative recombination continua are indicated above and below the spectra, 
respectively. 
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